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Abstract 

The conjunction of insights from structural biology, solution biochemistry, ge- 

PQ , netics and single molecule biophysics has provided a renewed impetus for the con- 

O ! struction of quantitative models of biological processes. One area that has been a 

r^ ' beneficiary of these experimental techniques is the study of viruses. In this paper 

we describe how the insights obtained from such experiments can be utilized to 

construct physical models of processes in the viral life cycle. We focus on dsDNA 

bacteriophages and show that the bending elasticity of DNA and its electrostatics 

^ ■ in solution can be combined to determine the forces experienced during packaging 

CN I and ejection of the viral genome. Furthermore, we quantitatively analyze the effect 

of fluid viscosity and capsid expansion on the forces experienced during packaging. 

\^ • Finally, we present a model for DNA ejection from bacteriophages based on the 

Cp , hypothesis that the energy stored in the tightly packed genome within the capsid 

Q I leads to its forceful ejection. The predictions of our model can be tested through 

"^ ■ experiments in vitro where DNA ejection is inhibited by the application of external 

• 1— I . osmotic pressure. 

> 

j^ . 1 Introduction 



Bacteriophage have served as a historical centerpiece in the develop ment of 

molecular biology. For example, the classic Hershey-Chase experiment ( Hershev and Chasd . 
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Fig. 1. Life-cycle of a bacterial virus. The ejectio n of the genome into the host cell 



happens within a m inute for phage like A and T4 (jNovick and Baldeschwielen . 
Letellier et al.l . 1200 4). The eclipse period (time between the viral adsorption a nd the 
first appearance of the progeny) is about 10-15 minutes ( Endv et al. . 1997^ . The 



packaging of the genome into a single capsid takes about 5 minutes (ISmith et al. 
200 L l. Lysis of the bacterial cell is completed in less than an hour (jEndv et al 



19971 '! 



19521 : lEchola . 120011 1 that establised nucleic acid to be the carrier of the genetic 



blueprint was performed using bacteriophage T2. The biology of bacterio- 
phage A provided a fertile ground for the development of the understanding 



of gene regulation (jPtashnd . l2004f ). while the study of virulent T phages (Tl 



to T7) paved the way for many other advances such as the definition of a 
gene, the discove ry of mRNA, and elucidation of the triplet code by genetic 
analysis (Echols, |2001). One of the central theses of the present paper is that 
phage are similarly poised to serve as compelling model systems for quanti- 
tative analyses of biological systems. Indeed, each of the stages of the viral 
life-cycle (see Figure 1) can be subjected to physical analysis. 

As shown in fig. 1 a typical phage life-cycle consists of: adsorption, ejection, 
genome replication, protein synthesis, self-assembly of capsid proteins, genome 
packaging inside the capsid and lysis of the bacterial cell. A wealth of knowl- 
edge about various aspects of these processes has been garnered over the last 
century with the main focus being on replication and protein-synthesis. How- 
ever, recent developments in the fields of x-ray crystallography, cryoelecton 
microscopy, spectrosc opy etc., have helped reveal the structu ral properties of 
the viral components ( Baker et al.l . llQQQt iLeiman et al.L 120031) inv olved in eiec 



tion, assembl y, and packaging , while recent experiments on 029 (ISmith et al. 



20011 ) and A ( Evilevitch et al.l . 120031 ) have helped quantify the forces involved 



in the packaging and ejection processes, respectively. In this paper we bring 



together these experimental insights to formulate quantitative models for the 
packaging and ejection processes. The model is based on what we know about 
the structural "parts list" of a phage: the shape, size, and strength of the cap- 
sid, the mechanical and electrostatic properties of DNA, and high-resolution 
images that reveal the structure of assembled phage particles. Our goal is to 
create a detailed picture of the forces implicated in DNA packaging and ejec- 
tion and, most importantly, make quantitative predictions that can be tested 
experimentally. 

The paper is organized as follows: in Section 2 we examine dsDNA bacterio- 
phage with the aim of assembling the relevant insights needed to formulate 
quantitative models of packing and ejection. In Section 3 we develop the model 
by examining the DNA packaging process in detail, and in Section 4 we show 
that it can be used to explain the DNA ejection process as well. The final 
section takes stock of a range of quantitative predictions that can be made on 
the basis of the model and suggests new experiments. 



2 Physical processes in the bacteriophage life cycle 



Since the early measurements of virus sizes many experiments involving viruses 
have been of a quantitative character. The emergence of quantitive insights 
into viruses has come from many quarters including electron microscopy. X-ray 
crystallography, single molecule biophysics and a large repertoire of classical 
methods in biochemistry and genetics. For many viruses, these techniques have 
given us a detailed picture of the structure and function of the entire viral parts 
list. Many viral genomes have been se quenced, and the structure of many o f 
the proteins coded for in the genome (JBaker et al.L llQQQt iReddv et al.L 120011 ) 



have been solved. Though it is impossible to discuss all of these advances 
here, we review some of the experimental insights that have informed our 
model building efforts. 



2.1 Experimental Background 



DNA is highly compressed inside bacteriophage capsids, and the resulting 
forces have important effects on the phage life cycle, as revealed in several 
experiments. One early experiment which shed light on the possible role of 



the fo rces associated with packaged DNA is that of lEarnshaw and Harrison 



(|1977|), who characterized the tight packaging of DNA in viral ca psids by the 



distance dg between strands (around 2.8 nm in full capsids), and bv lFeiss et al. 



(|1977|), wh o identified limits on the amount of DNA that can be packaged into 



a A capsid. iFeiss et al.l found an upper limit barely above the wild-type genome 



length, an d sugg ested that adding more DNA makes the capsid unstable. 
Rau et al.l (|l984r i made measurements on large volumes of non- viral DNA 
that showed that dg values in the range of 2.5 nm to 3.0 nm correspond to a 
pressure of several tens of atmospheres. 



Since the experiments of iFeiss et al.l and lEarnshaw and HarrisonL there have 
been a variety of experiments on DNA packing, several of which stand out 
because of th eir relevance to quan tifying the internal buildup of force dur- 
ing packaging. IShibata et al.l (|l987l ) measured the rate of packaging for phage 
T3, under various temperatures and chemical conditions. Through these ex- 
periments, they determined that the packaging process was reversible — one 
third of th e genome was ejected back into solution upon early interruption of 
packaging. ISmith et al.l (|200l[ ) reinforced the idea that a strong force builds 
up during packaging with real-time single-phage packaging experiments. They 
measured the rate of DNA packaging while subjecting the DNA to various 
resisting forces, quantifying the forces imposed by the packaging motor and 
the force resisting further packaging as a result of the confined DNA. Taken 
together, these two experiments give us a picture of DNA packaging in which 
a strong portal motor consumes ATP and in so doing pushes the DNA into 
the capsid against an ever- increasing resistive force. 



It is believed that the tightly wound DNA stores large energies resulting in 
high forces, which in some cases could aid DNA ejection into the host cell. Re- 



cently , an experiment to test this hypothesis was conducted bv lEvilevitch et al. 



(20031), who coerced A into ejecting its genome into a solution containing 
polyethylene glycol (PEG) to create an external osmotic pressure. They found 
that various osmotic pressures of several tens of atmospheres could halt the 
ejection process resulting in fractional genome ejection. The fractional ejection 
reflects a balance of forces between the inside and outside of the capsid. From 
this experiment it is concluded that forces are still present during ejection at 
the same high levels as were observed during packaging and in static capsids. 



There are a variety of impressive in vitro experiments which demonstrate 
pressure- driven ele ction of the phage genome ( Novick and Baldeschwieleii . ll988l : 
Bohm et al.L 120011 ). For the in vivo case, ejection driven by internal force is only 
one of several mechanisms that have been hypothesized to participate in trans- 
ferring the ge nome of bacterio phage into the host cell. Another mechanism is 
suggested by JMolineuxl (J200l|) who, on the basis of a wealth of experimental 
evidence, argues that the DNA of phage T7 is assisted into the cell by DNA 
binding proteins. It is likely that different bacteriophage use a combination of 
these two methods for ejecting their DNA into the host cell. 



Finally, recent cryo-electron microscopy (cryo-EM) and X-ray crystallography 
studies of bacteriophage have revealed their detailed internal structure and 
particularly, the ordered state of the packaged DNA. ICerritelli et al.l ()l997l ) 



verified the tight packing measured by lEarnshaw and HarrisonI and showed 
that the DNA is apparently organized into circular rings within capsids. Other 
structural experiment s have revealed t he structure of components involved 



duri ng ejection in T7 ( Kanamaru et all 12002) , the assembly of b acteriophage 



(^29(lTao et al.l . ll998l ). and structure of the packaging motor in (f)2 9 ( Simpson et al. . 



20001 ) ■ This structural information complements the single molecule measure- 



ments and will guide us in the construction of a quantitative model of the 
packaging and ejection processes. 



2.2 Orders of Magnitude in Bacteriophage Biophysics 



In the previous section we provided background on some of the experimental 
advances on bacteriophage which quantify the packaging and ejection pro- 
cesses. It is the aim of the remainder of this paper to discuss the implications 
of these experiments in a more quantitative way and to make predictions 
about the phage life cycle that can be tested experimentally. Before describ- 
ing our models in precise terms, we first perform estimates of the orders of 
magnitude of relevant physical properties involved in phage biophysics. Bac- 
teriophage r ange in size from a few tens of nanometers to several hundred 
nanometers ( Baker et al.L Il999| ) . The capsids of most are regular icosahedral 



structures a few tens of nanometers in size. Table 1 gives an idea of the typ- 
ical sizes of bacteriophage, as well as some animal viruses for comparison. 
These small containers house a genome which is several tens of microns long, 
a feat that demands extremely efficient utilization of space. In fact, a useful 
dimensionless quantity for characterizing the packaging efficiency is 

^genome r.. 

Ppack = -F, ' UJ 

^ ^capsid 



where ^genome is the volume of the genetic material and f^capsid ^^ ^^^ volume 
of the capsid. For double-stranded DNA bacteriophage, this result may be 
rewritten simply in terms of the number of base pairs in the phage DNA, 
A^]-,p, using the approximation that DNA is a cylinder of radius Inm and length 
0.34nm per base pair. Using these approximations, Ppack "^^^ ^^ rewritten as 

0-34^^bp 



Ppack Q 

^ 'capsid 



making the calculation of Ppack straightforward. Note that in this formula, 
^capsid i^ computed in units of nm^. 



Virus type 


Host type 


genome length (kbp) 


Diameter (nm) 


^pack 


Bacteriophage T7 


Bacteria 


40 


55 


0.490 


Bacteriophage (p29 ^ 


Bacteria 


19.4 


47 


0.459 


Bacteriophage T4 


Bacteria 


169 


92 


0.443 


Bacteriophage A^ 


Bacteria 


48.5 


63 


0.419 


Bacteriophage P22 


Bacteria 


41.7 


63 


0.319 


Herpes Simplex Virus 1 


Human 


152 


125 


0.159 


Human Adenovirus C 


Human 


36 


80 


0.143 


Smallpox Virus 1 ^ 


Human 


186 


220 


0.036 


Polyoma Virus SV40 


Human 


5.3 


~50 


0.083 


Mimivirus '^ 


Amoeba 


~800 


~400 


0.026 


Papillomavirus BPVl 


Animal 


7.9 


60 


0.070 



Table 1 

Packaged volume fractions of some bacteriopha ge and eukaryotic v iruses. We have 

used the outer dimensions of the capsid from iBaker et al.l ((l222l) in the calcula- 



tion of /Opapjj since these are more readily available than the dimensions of the 
empty space inside the capsid. The genome le ngths are given, for most viruses, in 
( National Center For Biotechnologv Informati on. 2004). The DNA in bacteriophage 
is seen to be significantly more tightly packed than the other viruses, revealing the 
geometric origin of the large packing forces associated with bacteriophage. 



For the purposes of examining the significance of this parameter, Table 1 
shows Ppack ^^^ ^ number of different viruses. A trend that is evident in the 
table is that viruses that infect bacteria are more tightly packed than the 
viruses that infect eukaryotic cells. A likely reason for this difference in degree 
of genome compaction is the difference in infection strategies employed by 
the two types of viruses. While eukaryotic viruses are brought into a host 
cell through processes in which both the genetic material and the capsid are 
taken into the infected cell, bacteriophage typically attach to the outside of 
the host and eject their DNA into the cytoplasm through a small channel. In 
order to trai isport their DNA quickly into the host, which itself is pressurized 
at ~ 3atm ( Neidhardtl . 119961 ). bacteriophage may power the ejection with a 
large internal pressure. However, as mentioned earlier, there is experimental 
evidence in the case of phage T7 that DNA binding proteins play a role in DNA 



(jTao et al.l ■ ll998^ Since the (/>29 capsids are aspherical, we use an average diameter 
th at gives the correct volume. 
Dokland and MurialdoL Il993[) 



World Health Organization! . 



200J) Since the smallpox particles are aspherical. 



we use an average diam eter that gives the correct volume. 

^ ( La Scola et al.l . 120031 ): this is the largest virus currently known. 



transport. These experiments raise doubts about the possibihty o f finding a 
single mechanism responsible for ejection from all types of phage (JMolineux, 

2nni[ ). 



The bacteriophage life-cycle is a dynamic process and it is of interest not 
only to consider the geometric parameters associated with viral DNA, but to 
attend to the temporal scales that are involved as well. The first step in the 
cycle is the adsorption of the phage onto the host cell. The frequency of this 
event depends on the abundance of available phage particles and their hosts. 
Since about 50 to 300 new phage particles are released by a single infection, 
the destruction of the host ce l l in a culture proceeds exponentially quickly 
once a cell is infected f Yound . 119921 1. After the phage has attached itself to 
the host its genome is generally r eleased on time scales ranging from seconds 
to minutes (|Letellier et al.L l2004 ). Probably concomitantly the transcription 
and translocation machinery of the host cell is hijacked and the production 
of phage proteins and factors required for replication of its genome begins. 
The time between the adsorption of the phage and t he appearance of th e first 
progeny capsids is usually on the order of minutes ( Flint et al.L l2000| ). This 
period is known as the 'eclipse period' and it is the time required to build 
up the concentration of the phage proteins to a level high enough to initiate 
self-assembly of the capsid, tail and motor proteins that constitute a mature 
phage particle. Self-assembly is a highly concentration dependent process, but 
once it starts it proceeds rapidly to completion in a few seconds. 



The second step, th e packaging process, is completed in about 5-6 minutes 
fjSmith et al.l . l200lh . The packaging rate is on the order of 100 bp/s in the 
initial stages but it slows down as more of the genome is confined inside the 
capsid. That is, the rate of packaging depends on the force opposing the motor 
as it packages the genome. This internal force grows as the amount of genome 
packaged increases, and the magnitude of the force depends on the solvent 
conditions. In fact, some biologically important multivalent ions, such as sper- 
midine, caus e spontaneous DNA con densation resulting in much smaller pack- 
aging forces (JEvilevitch et al.Ll2004 ). This will be demonstrated more clearly 
later in this paper. After the progeny phage have been completely assembled, 
enzymes lyse the host cell and a new generation of phage are r eleased. The 
proce ss of adsorption to lysis is completed in less than an hour (|Flint et al. . 
2000l) . 



The objective of this section was to highlight some of the opportunities for 
quantitative analysis provided by processes in the viral life-cycle. Despite the 
wide range of interesting physical process in the viral life-cycle, the remainder 
of this paper focusses on the physical forces associated with packaged dsDNA 
and the implications of these forces for the packing and ejection processes. 



3 The DNA Packaging Process 



The model we invoke to examine the energetics of packaged DNA is predicated 
upon two key physical effects: i) the elastic cost to bend DNA so that it will 
fit in the viral capsid and ii) the interaction energy which results from the 
proximity of nearby strands of DNA and which derives from charges on both 
the DNA and in the surrounding solution. 

O ur paper is very much i nspire d by previous theoretical work: The work 
of iRiemer and Bloomfieldl ( 1978^ laid the foundation for subsequent efforts 
on the energetics of packaged DNA by systematically ex amining each of the 



poss ible contrib u tions to the overall free energy budget. JKindt et al.l (|200l[ ) 
and Tzlil et al.l (120031) estimated the forces in packaging of the phage A. 
Arsuaga et al.l ( 20021 1 computed conformations of DNA in phage P4 using 
molecular mechanics models, observing that the conformation of DNA within 
a capsid depends on its volume, in agreement with our sugges tion that the 
parameter Ppa^k characterizes the extent of packing in phage. lOdiild (J2003|) 
has discussed several issues in bacteriophage packaging including the impor- 
tant problem of deriving an expression for the electrostatic interactions be- 
tween DNA strands starting from the Poisson-Boltzmann equation with the 
added com plication of high density of packing and po ssibly non-hexagonal ar- 
rangement, kindt et all (|200lMPurohit et al.1 (|2008bl ) and iTzhl et al.1 (|2003l) 
do not derive these expressions a priori; r ather they use results obtained 
by iRau et al.l ( 1984 ) , IParsegian et al.l ( 19861 1 , and iRau and ParsegianI ( 19921 ) 
from osmotic pressure experiments on DNA to characterize these interactions. 
We have already underlined the importance of electrostatic interactions in 
solutions as a c ritica l factor in the bacteriophage life cycle. More recently, 
Qdiik and Slokl ( 20031 ) analyze the possibility of a non- uniform density of DNA 
insid e the capsid and compute the size of regions within the capsid that are 
void. iMarenduzzo and Michelettil (20031) study the effects of the finite thick- 
ness of DNA on the forces experienced during packaging within a viral capsid. 
The aim of this paper is to take models like those described above and to sys- 
tematically examine trends from one virus to another as well as for mutants 
within a given phage type. 



3. 1 Structural Models for the Packaged DNA 



In order to construct a quantitative model of DNA packaging and ejection, 
it is necessary to specify the arrangement of the packaged DNA. One of the 
earliest successful attempts at determining the stru cture of DNA packaged 
i n a p hage capsid was an X-ray diffraction study bv JEarnshaw and Harrison 
()l977r ) on P22 and some mutants of A. Their studies revealed both long- range 



and short-range order of the encapsidated DNA. The evidence of long-range 
order conies from the observation that the diffraction is modulated by a se- 
ries of ripples indicating that the phage head is uniformly filled. Short-range 
order was indic ated by a strong peak corresponding to a 25A spacing. In 
an earlier study, iRichards et al.l ( 19731 ) used electron microscopy to visualize 
gently disrupted phage particles and found that the DNA close to the cap- 
sid boundary has a circumferential orientation. These investigations together 
suggested a coaxial spool like arrangement of the DNA inside the c apsid . 



This model wa s further corroborated by experiments of JBoov et al.l (|l99l[ ) 



on HSV-1 and ICerritelli et al.l (|l997|) who used cryo-electron microscopy to 
obtain three-dimensional visualizations of the packaged DNA in T7 phage. 
They consistently found a concentric ring-like geometry in both wild-type and 
mutant versions of T7 with a spacing of around 25A in nearly-filled phage 
heads. However, none of these experiments shed much light on the geometry 
of packaged DNA in the very early stages of packing. 



A well ordered structure of the encapsidated DNA is also suggested by other 
observations on DNA condensation in the presence of polyvalent cations or 
other condensing agents such as methanol, ethanol or dilute so lutions of sper- 



midin e, PEG and other low molecular weight polymers (see iGelbart et al. . 
2OOOI for an extensive overview). The strands in the condensate are known to 
have local hexagonal coordination, another piece of evid ence in favor of the 
kind of short-range order describ ed above. Sim ulations bv JKindt et al.l (|200l[ ) 
and theoretical investigations bv lOdiikI ( 19981 ) point in this direction as well. 
Evidence in support of this hypothesis is also provided by the simulations of 
Arsuaga et al.l ( 20021 ) who show that the coaxial spool (or inverse spool) is an 
energetically favorable configuration, particularly when the density of packing 
is high as in bacteriophage. In light of these arguments we assume that the 
DNA is arranged in columns of concentric hoops starting from the inner wall 
of the capsid. Each of the hoops is surrounded by six similar hoops (hexagonal 
arrangement) except those at the innermost column and those touching the 
surface of the capsid. Though we have assumed a limited class of geometries 
for the packaged DNA, namely spool models, we note t hat such model s are not 
necessarily the lowest possible free energy structures (|Klug and Ortiz . 2003) . 
Further, we remain unclear as to the precise dynamic pathways that would or- 
chestrate such structures during packing and ejection, though simulations are 



consistent with the dynamical development of su ch structures ()Arsuaga et al. 
2OO2I: iKindt et al.L l200l[ iLaMaraue eFID . HooJ 



3.2 Modeling the Free Energy of Packed DNA 



With a specific arrangement of DNA in hand, we can now compute the free 
energy required to package that DNA. Our efforts to write an expression 



for the free energy of the DNA packed inside a phage capsid are inspired 
by the experimental insights described above concerning the configuration 
of t he encapsidated DNA and it s behavior in ionic environments. We fol- 
low (iRiemer and Bloomfieldl . ll978l ). (|Odiild . ll998l ) and (|Kindt et alibOOll ) and 
break up the energetics of the encapsidated DNA into an elastic term and a 
DNA-DNA interaction term: 



G^tot(4, L) - Gbend + ^int • 



(3) 



We begin with the first term, by asserting that the bending energy in an elastic 
fragment of length L is given by 



G{L) 



ipksT 



dx 

W) 



(4) 



where ^p is the persistence length of DNA and R{s) is the radius of curvature 
at the position with arc length parameter s. This simply reflects the energetic 
cost G]2gj^(j = ^pkBTl/2R^ to bend an elastic rod of length / into a circular 
arc of radius of curvature R. 



The expression for the bending energy is considerably si mplified by a structural 
i nsigh t provided by t he exp erim ents described earlier o f JEarnshaw and Harrison 
( 1977t ). lBoovet all (|l99l[ ). and lCerriteUi et al.1 (|l997t) which showed that the 



DNA is arranged in the capsid in a series of circular hoops starting close to the 
surface of the capsid and winding inwards in concentric helices. Therefore, we 
specialize this expression to a hoop of radius R of length / = 2tiR and deduce 
that the bending energy in the hoop is Gj^QQp = n^pkBT / R. We ignore the 
pitch of the he lix and think of the en capsidated DNA as a series of hoops of 
different radii ( Purohit et al.l . I2003H ) . This leads to the following expression 
for the total bending energy: 



G'bend(^) = <P^B^II 



N{R, 

Ri 



(5) 



where N{Ri) is the number of hoops of radius Ri in the capsid. Note that in 
neglecting the helical pitch we have also assumed that the successive layers 
of hoops are parallel to each other. A more realistic model would acknowl- 
edge that strands in successive layers form a criss-cross pattern. We drop 
this complication in favor of the simple model used above that captures the 
essent ial physics. The p ersistence length of DNA depends on solvent condi- 
tions ( Smith et al.L 119921 1 and also on the sequence of base-pairs f Bednar et al. 
19951 ). However, ^p = 50nm is an appropriate nuni ber for the solvent conditions 



in the packaging experiment of lSmith et al.l ()2n0lh and the ejection experiment 



10 



of lEvilevitch et al.l ( 2003^ . We note here that the size of many bacteriophage 
capsids is on the order of a few tens of nanometers which is comparable to 
the the persistence length of double stranded DNA itself. This hints that 
the bending energy indeed constitutes a significant por tion of the free energy 
of th e encapsidated DNA, though recent experiments (|Cloutier and Widom . 
2004^ reveal that bending DNA for small radius of curvature may not be as 
costly as suggested by the linear elastic model used here. Note that we have 
neglected twist and w rithe as contributors t o the free energy y e ry much in 



hue w ith th e work of lArsuaga et al.1 JiooJ, iKindt et al.1 (|2nnil) . iTzlil et al 



f)2nn3t) . and lOdiikI (|2nn3l ). The twisting modulus of DNA is higher than its 



bending modulus and the twist would be relaxed if the ends of the DNA are 
free to rotate. However, there is no conclusive experimental data validating 
this claim and more work is required on this poi nt. There is a possibi lity that 
the DNA packaging motor is a rotational device ( Simpson et al.l . 120001 ) and in- 
vestigations elaborating its mode of operation are expected to shed more light 
on this problem. Given this state of our knowledge we neglect the twisting 
energy of the DNA and simplify the expression for the bending energy by con- 
verting the discrete sum in eqn. (5) to an integral according to the prescription 
Sj = -Jj^ / dR. In particular, we replace the sum of eqn. (5) with 



<^bend(-^) 



2'K^pkB T 
VSds 



R. 



out 



N{R') 
R' 



dR\ 



(6) 



where -Rout is the radius of the inner surface of the capsid and the length L 
of packaged DNA is given by 



R. 



y/3ds 



out 



R'N{R')dR'. 



(7) 



R 



The factor of ^ds appears because it is observed (jCerritelli et al.L Il997f ) that 
the DNA strands are arranged in a lattice with local hexagonal coordination 
and spacing dg. In other words, each strand has six nearest neighbors except 
those near the surface of the capsid and the innermost cylindrical space where 
there are three nearest neighbors on average. We expect that the interaction 
of the DNA with the pro teins of the capsid would give rise to surface energy 
terms ( Tzlil et al.l . 12003^ in the expression for the total free energy, but we 
neglect these terms in our analysis. 



Though not written explicitly in eqn. 3, entropy is another possible contributor 
to the overall free energy budget of the packaged DNA. We follow Riemer and 
Bloomfield (1978) in showing that this is small in comparison to the other 
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contributions to the free energy arising from the elasticity of the DNA and its 
interactions with itself. To see this we consider the extreme case of free DNA 
that can explore a large set of configurations in solution. This allows us to 
obtain an upper bound for the entropic contribution to the free energy, since 
we remove the constraint imposed by the tight packaging within the capsid. A 
fragment of DNA of length L in a solvent can be modeled as a random walk on 
a three dimensional lattice with lattice spacing 2^p,where ^p is the persistence 
length iDoil ( 19961 ). The total number of links, N, in the chain equals L/2^p 
and the total number of configurations of the chain on the lattice is z^, where 
z is the co-ordination number of the chain (6 in this case). Consequently, the 
free-energy contribution due to the entropy is given by 



G^ent 



LkT 



2ep 



ln6 



LkT 

IT' 



A similar expression for the entropic contribution can also obtained using 
Flory-Huggins type arguments following JRiemer and Bloomfieldl ( 19781 ). Since 
C,p = 50nm, for bacteriophage 029 with L = 6600nm, Gent ~ ISO/cbT. This is 
two orders of mag nitude smaller than the total energy of 2 x lO^ksT obtained 
experimentally bv lSmith et al.l (|200H ) for bacteriophage 029. We therefore ne- 
glect the entropy of the DNA in the rest of our calculations in this work. 



We next turn to the DNA-DNA interaction energy G^^^. DNA has a backbone 
that is highly negatively charged. As a result there are large energetic costs 
associated with bringing DNA fragments close together in solution. These in- 
teractions can in turn be screened by counterions which fill the space between 
DNA strands. From a theoretical perspective, it is natural to treat these inter- 
actions by applying the Poisson-Boltzmann approximation, in which charge is 
smoothly distributed acco rding to a Bo ltzmann distribution consistent with 
the potential it produces (jNelsonl . 120031 ). In this approximation, the free en- 
ergy is calculated from a sum of entropy (assuming the ions are locally an 
ideal gas) and electrostatic energy. Water is modeled as a continuous dielec- 
tric. In the limit of extremely close packing, ideal-gas pressure dominates over 
the relatively constant electrostatic e nergy, and the co rrect pressure is pre 



dieted. However, the experiments of iRau et al.l (|198J) and iParsegian et al. 



( 19861 ) show that this approximation does not give the correct dependence 
of the pressure on interstrand spacing as shown in Figure 2. The mismatch 
between theory and experiment is likely due to effects ignored by the Poisson- 
Boltzmann treatment, such as the discreteness of ions and of water molecules. 
However, detailed Monte-Carlo simulations in which the ions are treated as 



discrete objects but water remains continuous ()Lvubartsev and Nordenskiold , 
19951 ) work remarkably well. Though our treatment of the interaction energy, 
^int' ^ill t)e based upon the experimental data of JRau et al.l (|l984[ ). calcu- 
lations like those of iLvubartsev and Nordenskioldl ( 19951 ) could be used to 
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Fig. 2. Pressure in a hexagonal lattice of DNA, according to experiment and Pois 



son-B oltzmann theory. Experimental data points are from the data of iRau et al. 
((l^SJ) for 25 mM and 5 niM Mg ^ concentrations. Our theoretical calculations 
follow from a discrete one-dimensional Poisson-Boltzmann solver, assuming cylin- 
drical symmetry. The free energy was calculated as a sum of the Shannon entropy 
of the ions and the electrostatic energy of ions and DNA, with the zero point for 
the potential set so that internal and external ionic concentrations were related by 
a Boltzmann factor. The theoretical predictions differ from the experimental points 
by a factor of ten, though the slopes are approximately correct and it is difficult to 
distinguish between the data for the two different concentrations. Also shown is a 
least-squares fit to the empirical datapoints, resulting in the parameters c = 0.30 nm 
and Fo = 1.2 x 10^ pN/nm^. 

construct a "first-principles" model of the interactions, at least for monova- 
lent and divalent counterions. We rely on empirical data for this study but 
suggest the use of simulated data if experimental data is not available. 



The experiments of mau et al.l ()l984f ) and mau and ParsegianI ()l992f ) provide 
an empirical formula that relates osmotic pressure p to strand spacing d, 
the range of two to four nanometers as 



m 



p{ds) =Foexp( — -] 



(9) 



The value of the pressure is dependent on both the ion concentration and its 
charge. When working with this formula, it is important to remember that a 
large change in Fq can be mostly compensated by a small change in c, since 
most data points are in a small region far from the p-axis. Hence, even though 
c is relatively constant, we state the value used for both c and Fq together to 
avoid confusion. Figure 2 shows the experimental data for a solution containing 
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5 and 25 inM MgCl2 at 298 K, in which measurements reveal c = 0.30 nm and 
Fq = 1.2 X 10^ pN/nm^. These values should be appropriate for use whenever 
Mg^"*" ions are the dominant species and have a concentration of 5-25 mM, 
conditions satisfied in solutions commonly used in ph age experiments, includ- 



ing, for example, SM buffer and t he TM buffer used i n lEvilevitch et al.l ( 2003^ 



Unfortunately, the buffer used in lSmith et al.l (|2001[ ). includes enough sodium 
to suggest that the effect of Na^ must be considered in addition to Mg^"*". The 
solutions in vivo are far more complex and Fq is difficult to determine. Even 
solutions in vitro contain different concentrations of ions and determining Fq 
for each one of them through experiments would be impractical. 



The data on the measured pressure in terms of ds can be used to deduce the 
functional form of energy stored in the electrost atic interactions. We d o not go 
through the details here but refer the reader to lPurohit et al.l ( 2003bl ) for the 
full calculation. The calculation rests on the assumption of a pair potential 
interaction among A^ parallel strands of length / each packed in a hexagonal 
array with a spacing dg. The total interaction energy of this arrangement is 



Gint(^5 ds) = VSFo{c^ + cds)L exp( -] 

c 



(10) 



where L = Nl is the total length of the strands. This is the expression for the 
interaction energy when the ions in the ambient solution are monovalent or 
divalent cations. In this regime the interaction between the strands is entirely 
repulsive and these interactions decay as the strands are moved farther apart. 
For trivalent and tetravalent ions the physics is quite different. DNA i n suc h 
a solution can condense into hexagonally packed tori ( Raspaud et al.L Il998t ). 
In this regime there is a preferred spacing do between the strands, and the 
measurements are well fit by, 



Gint(^' ^s) = V^FoL (c^ + cds) exp 



:^^^-^)-{c' + cdo)-l{dl-dl)lin) 



where c = 0.14nm and do = 2.8nm. For ds < do the interaction is strongly re- 
pulsive and for ds > do it is attractive. This expression is a good representation 
of the free energy of interactions between the DNA molecules for sp acings dg 
less than or equal to the preferred value do ( Rau and Parseeianl . Il992l ) . In viral 
packaging we encounter interaxial spacings in exactly this range and hence we 
will use this free energy to study the effects of repulsive-attractive interactions 
on encapsidated DNA. 



Given that we have now examined the separate contributions arising from 
DNA bending, entropy and interaction terms, we now write the free energy of 
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Fig. 3. Idealized geometries of viral capsids. 
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^^^^dR' + VSPoic" + cds)L exp(-^) . (12) 
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An analogous formula holds for the repulsive-attractive regime. Note that this 
expression reports the free energy of the inverse spool configurations when a 
length L has been packaged by relating, R to L via (7). We will now show 
that the spacing between the strands varies in a systematic way during the 
packing process reflecting the competition between bending and interaction 
terms. 



3.3 DNA Spacing in Packed Capsids 



Our model makes a concrete prediction for the free energy Gtot of packaged 
DNA in any phage and for a wide range of solution conditions. In order to 
find Gtot for a particular phage, we minimize eqn.(12) by varying dg, under 
the constraint that L given by (7) is equal to the length of DNA already 
packaged. The expressions for N{R') will differ for different capsid geometries. 
Most capsids are icosahedral and we idealize them as spheres. Some capsids, 
e.g., 029, have a "waist". We idealize them as cylinders with hemispherical 
caps. Eventually, we will find that the geometry does not affect the overall free 
energy of packing as long as the internal volume of the idealization is the same 
for each geometry, once again reflecting the importance of the parameter Ppack, 
introduced earlier. Before we specialize to particular geometries (see fig. 3)we 
observe that N{R') = ^^^ where z{R') is the height of a column of hoops of 
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DNA situated at radius B! . Using this fact and differentiating (7) with respect 
to (is, while holding L constant, gives us ^ = —^^^fk- Minimizing Gtot with 
respect to the interstrand spacing ds gives, 

VSFoexp -4/c) = — — ^ ^— j^- ^ . (13) 

R^dl dl j«out R>z{R')dR' 



This equation represents a balance between the bending energy terms and 
the interaction terms. Note that if the size of the capsid is fixed then longer 
lengths of packed DNA imply smaller radii of curvature for the hoops since 
the strands want to be as far away as possible from each other for the case in 
which the interaction between the adjacent strands is repulsive. On the other 
hand, smaller radii of curvature eventually lead to large bending energy costs 
resulting in a trade-off that is captured mathematically in eqn. (13). 

We now specialize this result to particular geometries (i.e. particular choices 



of z{R')). For a sphere, z{R') = 2JR'^^^ — R'"^, and the equation which deter- 
mines the optimal ds reads 



V3Foexp(-rfs/c) = H ^ — 



R'dl dl ^Rl^^-R' 



[- 



+ ^2u^ ^out-v/i?Lt--R^.X . . 



For a cylinder with hemispherical caps, z{R') = h + 2JR^^^ — i?'^, where h is 
the height of the waist portion, and ds is the solution of the following equation: 



V3Foexp(-Vc) = ^ 



R^dl 
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Figure 4 shows ds as a function of the fraction of the genome packed for 
five different phage. We assume that all of them are packaged under the 
same (repulsive) condition s, with Fq = 2.3 x lO^pN/nm^ and c = 0.27nm 



as m 



Purohit et al. 



(2003b) cor r espon ding to the best visual fit to the data in 



the experiments of 



Smith et al.l (|2001l ) which was performed in a solution con- 
taining 5mM MgCl2 and 50mM NaCl. Table 2 gives the details of the geometry 
of the phage used in the calculation. 
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10 20 30 40 50 60 70 80 90 100 
Percent packed 

Fig. 4. The spacing dg in five different phage under fuhy repulsive conditions with 
Fq = 2.3 X lO^pN/nm^ and c = 0.27nni. These val ues of Fn and c resu lt in the best 
visual fit to the data in the experiment on (/)29 bv lSmith et al.l ( 20011 ) (5mM MgCl2 
and 50niM NaCl). The graphs show dg monotonically decreasing as more of the 
genome is packaged. T7 is the most tightly packed while T4 is most loosely packed. 
A and HK97 show an almost identical history of dg vs. fraction of genome packed, 
since the two are closely related structurally. 



Phage type 


Model 
Geometry 


^out (nm) 


h (nm) 


Genome 
length (nm) 


^pack 


T4a^asakiet,al.,2iMl) 


Capped cylinder 


39.8 


29.0 


57424 


0.442 


T7t rCerritelh et al.. 1997) 


Sphere 


26.6 





13579 


0.541 


629 fTao et al.. 1998) 


Capped cylinder 


19.4 


12.0 


6584 


0.461 


HK97 fLata et al.. 20001 


Sphere 


27.2 





13509 


0.503 


A (Baker et al, 1999) 


Sphere 


29.0 





16491 


0.507 



Table 2 

Idealized geometries of bacteriophage. The radius and height are determined by us- 
ing the volume available to the DNA. They have been calculated from experimental 
data about the geometry of capsids from several sources (in parentheses above), 
'"oack ^'-'^ some phage in this table is higher than corresponding values in table 1 
since this table uses internal volumes whereas table 1 uses the outer dimensions that 
are more readily available, f The T7 phage is unusual in that a part of its cylindri- 
cal tail (radius 9.5nm, height 28.5nm) protrudes into the empty space within the 
spherical capsid. The space occupied by this tail is not available to the DNA and 
we exclude it to get an effective radius. 



As noted earlier, a sufficient concentration of polyvalent cations suffices to in- 
duce effective attraction between adjacent DNA strands. Under such repulsive- 
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Percent packed 



Fig. 5. The spacing dg under repulsive-attractive conditions. We use Fq = 0.5 
pN/nm^, do = 2.8nni and c = 0.14nm. This correspo nds to a s o lution 



containing 5mM CofNH 3)fiCb, O.IM NaCl, lOmM TrisCl (|Kindt et al.l . 12001 



iRau and Parsegianl . |l99j). The spacing remains at the preferred value of 2.8nm 
for most of the packaging process, except in the end when volumetric constraints 
lead to smaller spacings as a consequence of high energy costs for maintaining this 
ds. 



attractive conditions the left hand side of eqn.(13) changes and we get 



A/3Fo(exp( 



do- ds. 






dl jl^o'^^R'z{R')dR'' 



(16) 



Here too, the idea is to solve for the optimal dg at each packaged length 
L. The results with the repulsive-attractive potential for different phage are 
shown in Figure 5. W e note the similarity in the observed trends with the 
results of Kindt et al.l (J200l|) who also study the packaging process with a 
repulsive-attractive potential. In the presence of surface energy the encapsi- 
dated DNA assurn es a toroidal configuration in the early stages of packing 
( Kindt et al.L 120011) much as it would have in solution in the absence of a cap- 
sid f Raspaud et al.l . 119981 1. At the later stages of packing the inv erse spool is 



the more energetically favorable geometry and the calculations of iKindt et al. 



show a smooth transition between these configurations. However, we assume 
that the inverse spool is the optimal geometry throughout the packaging pro- 
cess since we do not have surface energy terms. Note also that the trends 
in figure 5 are quite different from those with the fully-repulsive conditions, 
figure 4. Most importantly, in the early stages of packing under repulsive- 
attractive con^ition^th^_DNA^^traj]^ tend to be at the preferred spacing 
of 2.8nm ( Kindt et al.Ll200ll : lTzlil et al.L 120031 ). Volumetric constraints at the 
later stages of packing result in tighter packing and ds decreases to values 
lower than the preferred 2.8nm. Under fully-repulsive interactions the inter- 
strand spacing decreases monotonically as more DNA is packaged. 



For many bacteriophage, even were they to be packaged under conditions in 
which there is an effective attraction between adjacent DNA segments the 
genome would not be fit into the capsid if the interaxial spacing is at the 
preferred value of 2.8nm. This would result in repulsive interactions between 
adjacent DNA segments in the terminal part of the packaging process. Inter- 
estingly, in these circumstances, ds can be estimated using strictly geometric 
arguments. In particular, we equate the total volume available in the capsid 
with the volume of the packaged DNA. In particular, this implies 

Vcap = ^dlL, (17) 



where L is the length of the packaged DNA. This in turn implies 



* ^ f-WTL- '^«' 



or more generally ds oc 1/v L. In deriving this expression we ha ve assumed that 
the volume of the cylindrical void (u noccupied by the DNA) (JArsuaga et al. . 



20021: iKindt et al.Ll2nni[lodiikl . ll998t ) in the middle of the capsid is negligible in 



comparison to the volume of the capsid. The predicted scaling of ds with DNA 
length L above can be compared to several different experiments as shown in 
Figure 6. In par ticular, the sp acing has been measured bot h in A mutants 
dEarn s haw and H arrisonl Il977l ) and in the T7 bacteriophage ( Cerritelli et al. . 



19971 ). As seen in the figure, the scaling suggested by the model appears to 
provide a satisfactory description of the measured trends. More generally, note 
that the spacing of the packaged DNA is one of the key points of contact 
between our theory and experiment. In particular, both fig. 4 and fig. 5 are 
predictive in that they suggest how the spacing ds varies from one phage to 
another for different solution conditions. 



3.4 Forces During DNA Packing 



The phage genome is subject to various resistive forces during the packaging 
process. As discussed earlier, these forces depend on the size and geometry of 
the capsid, solvent conditions and the size of the genome. In this section we 
calculate those forces using the expressions for the elastic and the interaction 
energy derived in the previous section. 
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L"^''^ (Linnm) 



9.5 



Fig. 6. Comparison of measured spacings dg wit h —^ scaling law . The circles cor- 
respond to interstra nd spacing in T7 obtained by Cerritelli et alJ ( 19971 ). The stars 
represent data from lEarnshaw and HarrisonI ( 19771 ) for mutants of A-phage. We fit 
straight lines to both these data sets to show that the spacing scales with the inverse 
square root of the packaged length in capsids that are nearly full. 

3.J^.1 Forces from DNA confinement 

Using the model described above, it is possible to evaluate the forces resisting 
packaging for different phage. To obtain the force we differentiate the total free 
energy, Gtot) with respect to the packaged length L. Since the dependence of 
the spacing ds on L is already known we substitute it into the expression for 
Gtot rendering it a function of L alone. We then carry out the differentiation 
and obtain the following expression for the force: 



F{ds{L),L) = -^ = V3Foexp{c' + cd,) + ^, 



(19) 



where we assume that the radius R and the spacing dg are known functions of 
L. In Figure 7, we show the result of a series of such calculations for different 
p hage, all done assuming ionic conditions such as those used in the experiments 
of lSmith et al.l (J200l|) . The value of Fq a nd c corres ponding to these conditions 
were determined by a fit to the data of ISmith et a l. (2001) for bacteriophage 
$29. In principle more accurate values of these parameters can be obtained by 
a least-squares fit. The key point of Figure 7 is to illustrate the trends across 
different phage with particular reference to the way in which the maximum 
packaging force scales with Pp^ck- This is shown more clearly in Figure 8. 
We observe that the maximum resistive force scales roughly linearly with 
the packing density Pp^^k ^-cross different phage of varied shapes. We have 
also plotted the maximum resistive force for 029 in Figure 9 for different salt 
concentrations of the ambient solution. The repulsive interactions between the 
DNA strands grow progressively larger as the concentration of the Na"*" ions 
in solution decreases. This results in larger packaging forces as can be seen in 
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Percent packed 

Fig. 7. Comparison of forces during DNA packing process for different phage under 
fully repulsive conditions. T7 requires the largest force to package since it is most 
densely packed. T4 is at the other end of the spectrum requiring the smallest force. 
The data above corresponds t o Fn = 2.3 x 1 0^ pN /nm^ and c = 0.27nm obtained 
by a visual fit to the data of ISmith et al.l (|2nnil ) , who conducted the packaging 
experiment for phage 029 in a solution containing 5mM MgCl2 and 50mM NaCl. 




0.44 046 048 0.5 0.52 
Packing density (Ppg^) 



0.54 056 



Fig. 8. Maximum resistive force in different phage under three different repulsive 
conditions, (a) corresponds to Fq = 1.1 x 10^ pN/nm^ and c = 0.27nm, (b) cor- 
responds to Fq = 2.3 X 10^ pN/nm^ and c = 0.27nm and (c) corresponds to 
Fq = 3.3 X 10^ pN/nm^ and c = 0.27nm. The forces increase with the packing 
density Pp^ck ^^^ ^^^° with increasing Fq. 

Figure 9. 

We have also obtained the force required for packaging under repulsive-attractive 
interactions. The results can be seen in Figure 10. The forces are considerably 
smaller when the strands are at the preferred spacing do = 2.8nm. As with 
the DNA spacing discussed earlier the results in Figures 7 and 8 are predictive 
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Fig.9 

Fig. 9. Maximum resistive force in 029 for different salt concentrations. The maxi- 
mum force increases as the salt concentration decreases since the DNA interstrand 
repulsion becomes larger as the solution becomes more dilute. The values of Fq 
a nd c for the salt s olutions used in this figure were obtained from fits to the data 
of JRau et al.] (|l984l ^. 




10 20 30 40 50 60 70 80 90 100 
Percent packed 

Fig. 10. Comparison of forces during DNA packing process for different phage under 
repulsive-attractive conditions with Fq = 0.5 pN/nm^, do = 2.8nm and c = 0.14nm. 
This co rresponds to a solution containing 5mM Cof NH-^lfiCb, O.IM NaCl, lOmM 
TrisCl ( Kindt et al.l . l200ll : iRau and Parseeianl . Il992h . The trends seen here are no 
different from those in the fully repulsive conditions - T7 requires large forces and 
T4 requires small forces for packaging. The maximum force, however, is significantly 
smaller than that seen for fully-repulsive conditions. 

and suggest a wide range of new experiments using both different solution 
conditions and different phage. 

The history of force as a function of packaged length obtained above rested 
on the assumption that the integral approximation adopted in eqn.(6) is an 
accurate representation of the bending free energy. In principle the optimal 
spacing ds{L) and the resistive force F{L) can be obtained by minimizing the 
free energy without resorting to the integral approximation. This has been 
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Fig. 11. Force and inter axial spacing as functions of the amount of DNA packed 
in bacteriophage (i)2'd. The hexagons correspond to the experimental data of 



m bacteriopnage 0z\). ine nexagons correspond to tne experimental aata or 
Smith et al.l (200l|). The thick lines are results of the continuum model and the 



circles connected with the thin line are obtained from the discrete model. 



carried out in lPurohit et al.l (J2003a!) and the results are shown in Figure 11. 
There are several distinctive features to be noted in this figure. Foremost 
among them is that the curves for the history of resistive force and interaxial 
spacing as a function of the length of DNA packed are not monotonic, un- 
like the curves obtained from the integral (continuum) models. The discrete 
steps represent the addition of a new stack of hoops, during the packaging 
reaction. More importantly, discreteness of packing implies that at specific 
lengths packaged there will be sharp changes of the spacing d^ due to an in- 
crease in the number of hoops per layer. These events might be a mechanical 
signal of the packaging configuration. In particular, the observation of such 
steps in an experiment would be evidence in favor of our quasi-static picture 
of the dynamics of packaging. 



3.4-2 Forces due to viscous dissipation 

Thus far we have calculated the free energy stored in the compressed DNA 
within a bacteriophage capsid, which represents the total amount of work that 
the motor must do to package the entire genome. We have not yet considered 
irreversible work that may result from the high rate of packaging. An obvious 
source of irreversible work is viscous dissipation in the fluid. Her e we demon 



strate that, in fact, such forces are negligible, as speculated by ISmith et al. 

( iiih . 



We identify four sources of viscous dissipation during packaging: drag on (i) the 
capsid and (ii) the unpackaged genome as each is pulled through the fluid, (iii) 
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R 


Capsid radius 


30 nm 


tc 


Capsid thickness 


2 nm 


Rp 


Pore radius 


2 nm 


N 


Number of pores 


10 


X 


Sheath length 


50 nm 


D 


Inner sheath radius 


1.3 nm 


d 


DNA radius 


1.0 nm 


V 


DNA velocity 


100 bp/s « 30 nm/s 


Rg 


Unpackaged genome radius of gyration 


300 nm 



Table 3 

Physical numbers for DNA packing in 029, taken from lTao et al.l (|l99a l. 



viscous dissipation within the sheath (as might be important during ejection), 
and (iv) dissipation as fluid is extruded through capsid pores. We calculate 
an upper bound on each to provide an upper bound on fluid dissipation in 
general, and conclude that such forces are negligible. Note, however, that we 
assume water to behave as a continuum material throughout this analysis. 
Non-continuum effects, which could affect our conclusions, are not treated 
here. In what follows, we use values from Table 3. 



First, as DNA is pushed into the capsid, equal and opposite forces pull both the 
unpackaged genome and the capsid through the surrounding fluid. Although 
the (smaller) capsid certainly moves more quickly than the genome, as an 
upper bound, we assume each to move at the full translocation velocity V. 
Furthermore, an upper bound can be obtained using the Stokes drag F = 
GnfiRbV on a solid sphere of radius Rb bounding each. Using Rb ^ R for the 
capsid and Rb = Rg for the genome, we obtain 



-^ cap 



-4. 



2 X 10"^ pN and Fgen ~ 2 x lO^^pN 



(20) 



jiving power dissipation 

[/cap ~ 1.5 X 10~^A;bT/s and [/gen ^ 1-5 x IQ-^ksT/s. 



(21) 



Second, to estimate the dissipation within the sheath, we consider the fluid 
drag on a DNA molecule (modelled as a cylinder) of radius d moving at velocity 
V through a cylindrical sheath of inner radius D and length X, into or out of 
a viral capsid. The fluid between the DNA and the sheath is assumed to obey 
the Stokes equations subject to the no-slip boundary conditions. Packaging or 
ejecting DNA requires an equal volume of fluid to be expelled from or injected 
into the capsid, which can occur either through pores in the capsid, or back 
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through the sheath. Below we treat the fluid dissipation in both cases. 



This is a textbook problem in fluid mechanics iLandau and Lifshita ( 1987 ) 



giving a fluid velocity profile (with radial coordinate r) 



The pressure AP depends on the nature of the capsid. If the capsid is imper- 
meable to water, the volume of DNA ird'^V entering the capsid must exactly 
equal the fluid volume 27r J^ urdr leaving the capsid, which gives 

AP = 4^^X 

{{d^ + D^)ln{D/d) + d^-D^)' ^ ^ 



of order 10^ Pa. The total power dissipated comes from shear stress on the 
DNA, t/shear = ^TidXV dru\r=d and from p — V work, driving a volume flux 
of DNA against a pressure AP, f/pressure = T^d'^VAP. These two give a total 
dissipation 



which represents an upper bound for dissipation within the sheath. 

If, as we assume below, it is easier for the fluid to flow through capsid pores 
than through the sheath, the 'backflow' in (22), proportional to AP, disap- 
pears. In that case, only shear stresses occur and give a total dissipation 

t/^ath = 2^^!^^^^ ~ 3 X l^-'kBT/s. (25) 



Finally, we estimate the power dissipated for fluid extruded through capsid 
pores, rather than the sheath. These pores occur at N symmetry points on 
the capsid shell. A crude estimate for the dissipation through each is obtained 
by assuming Poiseuille flow through each pore. Entrance and exit effects con- 
tribute, at most, a term of comparable magnitude. The flow rate through A^ 
pores due to a pressure AP inside the capsid is 

NAPttRI , ^ 



Requiring the DNA volume entering the capsid ticP'V to equal the fluid flux 
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Qp out of the capsid determines AP to be 

AP ~ ,;^, , 27 



giving rise to a viscous dissipation 

U, = APnd'V ~ ^^^^^ ~ lO-^fc^T/s. (28) 



Since the dissipation through the pores is so much lower than that through 
the sheath, one expects most fluid extrusion to occur primarily through the 
pores. 

All of these effects are negligible in comparison with the power supplied by the 
motor, which consumes roughly ~ 10k bT per 2 base pairs, giving PVmotor ~ 
400fcBT/s. 

Lastly, we note that if viscous dissipation is the primary damping mechanism 
during ejection through the sheath, the DNA ejection velocity can be obtained 
using the above results. If fluid must flow through the sheath to replace the 
volume lost by ejecting DNA, the ejection rate Ve ~ 6/im/s per atmosphere of 
applied pressure; whereas if fluid can flow freely through capsid pores, higher 
ejection rates Vg ~ 300/im/s, per atmosphere of applied pressure. Under the 
first assumption, 10-lOOkbp genomes (3-30 fim) would be ejected in less than 
a second, and fifty times faster under the second. 



3.5 Capsid Mechanics 



In all the calculations above we have assumed the capsid to be rigid. The 
capsid is actually a deformable object and it undergoes significant changes in 
shape and size during the DNA packaging process. These inelastic changes 
to the capsid geometry occuring at the early stages of packing are referred 
to in liter ature under the general headi ng of capsid maturation or prohead 
expansion ( Blacki . llQSSt iLata et al.L 120001 ). The exact mechanism behind these 



changes are not comp letely unders tood and are probably not due to the forces 



exerted by the DNA (jBlackl . 1198 81. In this section we will be concerned with 
the later stages of packing after maturation in an effort to understand the 
interplay between the forces exerted by the packaged DNA and the elastic 
deformation of the mature capsid. In what follows we treat the capsid as an 
elastic shell to obtain estimates about the maximum internal pressure it can 
sustain. We also examine the effect of capsid elasticity on the packaging forces. 
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The mature capsid of bacteriophages such as A, Tl and HK97 is made up of 
copies of a few proteins arranged on an icosahedral shell. For the purposes of 
this analysis we model the capsid as an elastic sphere and assume that the 
protein sub-units making up the capsids undergo only small deformations so 
that the elastic energy stored in a capsid expanded to a radius -Rout is given 
by 

Ccapl-Rout) = '^(-^out - Ro) , (29) 



where Rq is the equilibrium radius of an empty capsid and k is a constant 
measuring the ca psid stiffness. An examp le of such an energy emerges from 
the simulations of lTama and Brooks! (J2002) who perturb the positions of atoms 
in the capsid of a plant virus (CCMV) and measure its energy as a function of 
radius. They fit a quadratic polynomial in -Rout ^^ ^^e energy obtained from 
their simulations and find 

K = 5.0 X lO^pN/nm. (30) 

The capsid can be described as a thin shell of radius -Rout- I^ ^^e presence of 
an internal pressure p, the free energy is minimized when 

- ^nRl^^p + 2K{Ront - Ro) = , (31) 

which implies 

27!-R'^^^p 
-Rout - -Ro = ■ (32) 

K 

Using eqn.(30) and parameters from a typical phage capsid, p = 40atm and 
-Rout = 30nm, we find 

^out - ^0 = 4.6 X 10^2 nm . (33) 

This change is negligible in comparison to -Rout implying that a phage capsid 
can be treated as a rigid shell. 

The analysis given above estimates the deformation of the capsid in response 
to the forces exerted by the packaged DNA. It is also of interest to deter- 
mine the maximum pressure that a capsid can sus tain particularly in vi ew of 
the osmotic shock experiments on these systems ( Cordova et al.l . 120031 ). We 



approach this problem by modeling the capsid as an assemblage of proteins 
interacting through weak forces such as van der Waal's forces and hydrogen 
bonds. We note that capsids have thin walls compared to their diameter. For 
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example, the capsid of 029 is ab out l.Snm thick while its linear dimensions 



are of the order of 40-50nm (see IXao et all Il998l for data on 029). We use 



these ideas in conjunction wit h a coarse-grained model for the cohesive ener- 
gies between protein subunits (jReddv et al.1 . 120011 ) and estimate the maximum 
pressure sustainable for a ca psid to be in excess of lOOatm. The details of the 



calculations can be found in lPurohit et al.1 ( 2003bl ) 



4 The DNA Ejection Process 



We saw in the previous sections that packaged bacteriophage capsids are pres- 
surized with pressures as high as 60atm. This has led to the speculation that 
the high pressure in the b acteriophage provides the drivi r ig force for DNA 
electi on into the host cell (jSmith et all 1200 It JKindt et all 1200 It JTzlil et al. , 
2003^ . In this section, we examine the feasibility of t his hypothesis. Specifica lly, 



we show that internal forces explain the results in lEvilevitch et al.l ( 2003^ 



on 



inhibition of DNA ejection from phage A and allows us to make predictions 
for bacteriophages with varying genome lengths. 



Any experiment in which we can control the amount of DNA ejected from a 
bacteriophage by the application of external pressure can help us understand 
wheth er internal forces drive ejection. In a recent experiment bv lEvilevitch et al 
( 2003^ . A bacteriophage were coerced into ejecting their DNA in vitro with 



the help of a protein called LaniB or maltoporin. This protein, found on the 
outer membrane of E. coli is the natural receptor for A. When the phage 
binds to this protein it ejects its cargo of DNA. The DNA was ejected into 
solutions of polyethylene glycol 8000 (PEG) of v arious concentrat ions, which 
applied known osmotic pressures to the capsid. lEvilevitch et al.l found that 
osmotic pressures of 20 atm were sufficient to prevent any DNA from leaving 
the capsid, whereas DNA was partially ejected at lower external pressures. 

We can use our model for the forces associated w ith the packaged DNA to 
analyze the experiments of lEvilevitch et al.l ( 20031 ). To that end we consider 
the system of bacteriophage, ejected DNA and the PEG solution at equilib- 
rium. The energetics of the DNA inside the phage capsid remains the same 
as discussed earlier in this paper. The extra feature we add is the free energy 
of the ejected DNA-PEG system. Since the persistence length of the DNA 
(50nm) is much larger than the persistence length of the PEG molecule (~ 
1 nm) we model the insertion of the DNA inside the PEG solution as be- 
ing equivalent to the insertion of a rigid cylindrical rod of radius R inside 
a solution whic h exerts osmotic pressur e Hq on the rod. This problem has 
been studied bv ICastelnovo et al.l (|2003 ^. who estimate the work of insertion 
of rigid rod (DNA) into a polymer (PEG 8000) solution as a combination of 
pressure-volume work, energy associated with creating new surfaces and the 
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entropic effects associated with polymer in solution. Here, however, we resort 
to a simple approximation where we retain only the term associated with the 
pressure- volume work. Hence the work expended to insert length Lq — L oi 
DNA of radius R into the PEG solution is given by 



w{Lo -L) = Ho(Lo - L)txR^. (34) 

The total free energy of the system is the sum of the free energy of the DNA 
inside the phage capsid and the work of insertion and is given by 

Gtoi{ds{L), L) + Ho(Lo - L)TiR\ (35) 



We already know that the free energy of the DNA inside the capsid depends 
on the parameters Fq and c which in turn depend on the i onic strength of 



the bu ffer used in the ejection reaction. The experiment by lEvilevitch et al 
( 2003^ involved a buffer of 10 mM MgS04. Thus we can use the empirical val- 
ues c = 0.30 nm and Fq = 1.2 x 10^ pN/nm^ for Mg^^ solutions as determined 
earlier. We will use these values for all the fits to the experimental data and 
for further predictions. 



In order to find find Lq — L, the ejected length, we need to minimize the free 
energy Giqi with respect to L. Differentiating eqn.(35) we get, at equilibrium. 



dGtot{ds{L),L) ^ _. 



dL 



- lioTiR = 0. (36) 



The first term in the above expression is merely the resisting force F{d{L), L) 
derived earlier. Consequently equation(36) becomes 



Fid,{L),L) = UonR^. (37) 

This equation is solved for L for several values of Hq and the results for the 
percentage of DNA ejected as a function of osmotic pressure are given in Fig- 
ure 12, for different total DNA lengths Lq. Experiments are currently underway 
to examine the extent to which the predictions described here are borne out 
experimentally. Figure 13 shows how the model applies to th e ejection behav- 



i or of other phage under the same solvent conditions as in lEvilevitch et al. 
(2Qq3). 
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Fig. 12. Fractional DNA ejection in A-ph age as a function of osmotic pressure cor 



responding to experimental conditions in lEvilevitch et al.l (|20O3) (lOniM MgS04V 



The i nset shows the best visual fit to the experimental data of lEvilevitch et al. 
(|2nn,'^ for A with genome size 41.5 kbp usin g parameters -Fn and c, and also a pa- 



rameter free curve obtained from the data of iRau et al.l (| 19841 1 . We see a good match 
between the experimental findings and the theoretical results. Also shown are the 
predictions for ejection behaviour of A-phage for genome sizes of 37.4 Kbp, 45 Kbp, 
48.5 Kbp and 51 Kbp under similar experimental conditions. 
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Fig. 13. DNA ejection as a function of osmotic pressure for various wild-type s pecies 



of bac teriophage for ionic conditions similar to the experiment of Evilevitc h et al. 
(J2nn3h (lOmM MgS04). The hues show the DNA ejection behavior for T4, HK97, 
(j)29, T7 and A. Osmotic pressures as high as 35-55atm are required to inhibit ejection 
in these bacteriophage. 

5 Conclusions 



This paper addresses physical processes in the viral life cycle through a quanti- 
tative framework based on insights from structural biology, single molecule bio- 
physics, electron microscopy and solu tion biochernistry. The mo dels were moti- 



vated by specific experiments on 029 (|Smith et al.L 120011 ) and A (|Evilevitch et al. 
120031 ) but their applicability extends to all dsDNA bacteriophages. In fact, we 
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use our models to predict important features of the packaging and ejection 
processes in phages other than 029 and A. 

The key predictions arising from the modehng efforts described here are: 

• Dependence of Forces and Spacing on Ionic Strength. As shown in figs. 4,5,7 
and 10, there is a strong dependence of both the spacing of the packaged 
DNA as well as the forces that build up due to packing on the ionic condi- 
tions during the packaging reaction. We suggest systematic experiments to 
explore these effects. 

• Dependence of Forces and Spacing on Phage Identity. As shown in figs. 7 
and 8, we find a systematic and strong dependence of the packing forces on 
the particular phage species of interest. In particular, had the experiment of 
Smith et al. been carried out in phage T7, we predict a maximum packing 
force in excess of lOOpN, rather than the 57pN found in 029. A simple 
parameter for developing intuition concerning the forces associated with 
different phage is ppack-, the ratio of the volume of the genome to the volume 
of the capsid. 

• Force Steps During Packaging. One of the weakest points of the analysis 
described in this paper is the uncertainties that attend the particular struc- 
tural arrangements of the DNA on the way to the fully packaged state. In 
particular, we have assumed a sequence of structural states which are all of 
the inverse spool form and one consequence of this structural picture which 
might be testable is the presence of steps in both the DNA spacing and 
forces as shown in fig. 9. 

• Dependence of ejection inhibition on genome length, virus type and solu- 
tion conditions. The beautiful experiments of lEvilevitch et al.l ( 2003^ pro- 



vide a direct window on the forces associated with the packaged DNA. 
Fig. 12 and 13 represent a wide range of parameter-free prediction for the 
fractional ejection inhibition that should be seen in such experiments. 

Finally, a ca utionary note. Ex periments have shown that some bacteriophage. 



such as T7 (jMolineuxl . 120011 1 may rely on a different mechanism for deliv- 
ering their genome into the host cell. This is a possibility worthy of further 
exploration, but we emphasize that viruses may use many different methods 
or combinations thereof to propagate themselves. 
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